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Calculated Energy Dissipation Distribution in Air by Fast 

Electrons From a Gun Source x 

John E. Crew 2 

(November 22, 1960) 

Results of calculations on the energy dissipation distribution for electrons from a point 
collimated (gun) source in an infinite air medium are presented. The calculation has been 
made for a monoenergetic source of 0.4 Mev electrons. The method of moments has been 
employed, fitting the two spatial variables separately. 



1. Introduction 

Calculations of the energy dissipation distributions 
r for fast electrons in infinite and homogeneous media 
have been reported for the simplest source geom- 
etries. 3,4 The problem of electron penetration is 
considered to be the following: given a source of 
electrons in a material, calculate the energy deposited 
by the electrons in a small spherical volume as a 
function of its position in the material. Spencer has- 
calculated the energy dissipation distributions for 
the plane perpendicular source and the point iso- 
tropic source, 3,4 using the moment fitting technique 
which lias proved successful in treating X-ray 
penetration. 5 His calculations were made for a 
variety of materials and for a range of source energies 
from 25 kev to 10 Mev. Agreement with available 
experimental data has been good. 

The present work is concerned with a source 
geometry involving two spatial variables, namely 
that of the point collimated (gun) source, and 
corresponds to the following experimental situation: 
A collimated beam of monoenergetic electrons 
originates at a point in the medium. As the electrons 
move about, they dissipate energy. The quantity 
we determine, J(r,a), is basically the energy dissi- 
pated in a small volume at a distance r from the 
source point, at an obliquity angle a with respect to 
the initial line of fire, and at an azimuthal angle </> 
relative to an arbitrary reference plane. 

2. Calculation of the Moments 

Following the notation of footnote (4), the elec- 
trons are assumed to be generated at energy E , 
with residual range r and stopping power (—dE/ 
dr) E o. If N electrons are produced by the source, 
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then the energy dissipated in a small volume dV=r 2 
drd<pd(cosa) at the position (r, <*,</>) is 



dE N(dE/dr) E0 7 , , 
, 2 -J(r,a). 



dV 



?o 



(1) 



Determination of the function J(r y a) is our objective. 
This is accomplished by fitting moments of the 
closely related function 



F(r,a)=2ic(r/ro) 2 J(r,a). 



(2) 



F(r,a) is proportional to the energy dissipated per 
unit volume in a narrow ring of the material, whose 
symmetry axis is the initial line of fire, at a position 
(r,a) . The moments 

F n>l =- ( r °(r/r ) n dr P d(cosa) P^cosa) F(r,a). (3) 
r Jo J-i 

Thus for total energy production by the source 
NE , the lowest order coefficient is 



^o,o=4 fj(r,a)dV~ 
'o J 



E n 



r (dE/dr) E o 



(4) 



The F n j may be represented as linear combinations 
of a set of basic coefficients 

2w\ dx x" d(cosa) P t (cosa) Ii (t,a,x) (5) 
obtained by solution of the differential equation 

— dt +cosa 'bx =z \ *' J ^( cosa ')S0,e) 

ll l0 (t,a',x)-l l0 (t,a,x)} + ( ^p- S(x)8(t-l) P^cosa). 

(6) 
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This equation describes penetration by electrons 
originating on a plane, with intensity in different 
directions given by a Legendre polynomial. The 
parameter t is the remaining fraction of the electrons' 
initial range, x is the distance from the source plane 
in units of the initial electron range, and the scatter- 
ing kernel S(t,9) has an energy dependence expressed 
by a factor [tit+a)]' 1 . (See footnote [4].) 

Using the stopping power representation of eq. (14) 
of footnote (4) together with eq. (12.19) of footnote 
(5), which by symmetry applies also to the case of 
isotropic detector and collimated source, we write 



n ' l ~ (2J+l)n!(J)! £^« J ».« 



(n+i-%), I 



(7) 



where the ^4/s are given in footnote (4). 

Spencer has kindly made available to us tabula- 
tions of In;o° obtained with a generalized version of 
the machine program described in footnote (4). We 
have calculated from these a triangular set of F n j 
values with l<n<\2—l, 0<Z<6. The spatial 
moments F nJ were used to compute a set of functions 
F t (r) by the function fitting procedure of footnote (4). 
These moments are listed in table 1. Computations 
have been carried to six digits throughout for the sake 
of internal consistency. 

Table 1. Spatial moments for a point collimated source of 0.4 
Mev electrons in air 

The entries written as N(m) are to be interpreted as iVX10 n 



n 


I 


F n ,l 


n 


I 


F n ,l 



1 
2 
2 
3 



1 


2 

1 


0.143977 (1) 
.512249 (0) 
.383407 (0) 
.176665 (0) 
.182270 (0) 


9 
9 
9 
9 
9 


1 
3 

5 

7 
9 


0.232118 (-1) 
.838418 (-2) 
.157837 (-2) 
.173430 (-3) 
.123226 (-4) 


3 

4 
4 
4 
5 


3 


2 
4 
1 


.558737 (-1) 
.156178 (0) 
.741548 (-1) 
.162271 (-1) 
.827478 (-1) 


10 

10 
10 
10 
10 



2 
4 
6 

8 


.217520 (-1) 
.116838 (-1) 
.309866 (-2) 
.457991 (-3) 
.414279 (-4) 


5 
5 
6 
6 
6 


3 
5 

2 
4 


.265519 (-1) 
.436366 (-2) 
.744594 (-1) 
.369759 (-1) 
.854666 (-2) 


10 
11 
11 
11 
11 


10 
1 
3 

5 

7 


.200647 (-5) 
.134557 (-1) 
.511519 (-2) 
.103639 (-2) 
.123527 (-3) 


6 

7 
7 
7 

7 


6 
1 
3 

5 

7 


.109575 (-2) 
.422286 (-1) 
.144003 (-1) 
.251286 (-2) 
.258849 (-3) 


11 
11 
12 
12 
12 


9 
11 

2 

4 


.937227 (-5) 
.539047 (-6) 
.127403 (-1) 
.706332 (-2) 
.199041 (-2) 


8 
8 
8 
8 
8 



2 
4 
6 

8 


.389695 (-1) 
.201798 (-1) 
.500374 (-2) 
.683193 (-3) 
.578889 (-4) 


12 
12 
12 
12 


6 
8 
10 
12 


.317900 (-3) 
.312289 (-4) 
.201874 (-5) 
.922263 (-7) 



F l (r) 



= 2irp F(r,a 



3. Construction of the Energy Dissipation 
Function 

The fitting of the moments F n j in r yields a set of 
functions 



The same asymptotic form for the trial function for 
Fi(r) was used as for the plane perpendicular source 
calculation in footnotes 3 and 4. The function 
Fiir) is normalized to unity at the origin, this 
boundary condition being assured by choosing the 
trial function 

i^(r)=(l-£) *exp [-A l r/(r -r)} 



=0 



for0< - <0, 
~ r ~ 

for->l. 



(9) 



The asymptotic constant A t is calculated from the 
two highest moments of F n>i , and the coefficients 
ctifii are computed by fitting moments according to 
the procedure of footnote (4). The parameter 7 was 
set equal to zero for 0</<3. The choice 7=1 gave 
a better fit for 1=4. The results of the fitting in r 
are shown in figure 1 . The results for 1=0 correspond 
to the case of a point isotropic source and are in 
close agreement with the previous calculation by 
Spencer. 4 

To obtain F(r,a) it is necessary to sum the Legendre 
series 

F(r,a) = ±,(l + }i)Fi(r)P t (cos a) (10) 

1=0 

using the set of Fi(r)'s known for 0<Z<4. How- 
ever, the series has very poor convergence unless 
terms with Z>>4 are included; hence it was necessary 
to perform an extrapolation to higher I values. 
By fitting the F x to a suitable form and making use 



) P; (cos a) d(cos a). (8) 




r/r , FRACTIONAL RANGE 

Figure 1. FirU jive partial moments in r versus radial 
distance from the source in units of the fractional range, 
r/r . 

Fi (r) is in units of stopping power of 0.4 Mev electrons in air per unit range 
per electron. 
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of the generating function for Legendre polynomials, 
it was possible to sum the series. 6 The following 
form was used over an appreciable portion of the 
electron range: 



^i=S aid+}i) n exp l-bS+H)] 



(ID 



for it an integer. If in the generating function for 
Legendre polynomials 



(l-2xz+ ^)- l/2 =^,z l Pi(x) 



1=0 



(12) 



the substitutions z=e a and z = cos a are made, it is 
easy to show that 

F( r ^)^=J2 a>t - ' avhT ( cosn b i~ cos a )~ 1 ' 2 - 



V2 bW 



(13) 



fl L. V. Spencer, private communication. 
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Figure 2. Energy dissipation function versus angle with 
respect to initial beam direction for various radial penetra- 
tions, r/ro. 

J (r,a) is in units of stopping power of 0.4 Mev electrons in air per sterradian 
per unit range per electron. 



This method was used for r<0.5r with n=l. For 
larger values of r the series converges more rapidly, 
so that a simpler method of extrapolating to higher 
Fi values was used. A plot of log [i^/(/+l/2) w ] versus 
I was made. A value of m yielding a straight line 
for the last three or four calculated Fi was chosen by 
trial and error and the higher values obtained by 
extrapolating this line. The series was then summed 
explicitly with only a moderate number of higher I 
terms required. From F(r,a) the desired energy 
dissipation function J{r,a) was obtained according 
to eq (2). A plot of the energy dissipation function 
versus a for different penetrations r/r is presented in 
figure 2. A plot of J(r,a) versus r/r for different 
angles a is given in figure 3. A polar diagram 
showing contours of constant energy dissipation is 
given in figure 4. 
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Figure 3. Energy dissipation function versus radial penetra- 
tration for various angles with respect to initial beam direc- 
tion. 
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Figure 4. Contours of constant energy dissipation in air for 
0.4 Mev electrons. 

Contours of 2ttJ (r,a) are shown, with the radial distances giving the penetra- 
tion! n cm of air at 76 cm of Hg and 25 °C. 



(Paper 65A2-93) 



4. Discussion 

At present there are no experimental data with 
which to compare our results. Results of calcula- 
tions for the one-dimensional problem have yielded 
results in good agreement with experiment. Calcu- 
lations with slightly different trial functions generally 
reproduce results within about 3 percent. Since 
we are fitting twice, the errors can be expected to 
multiply, and our results cannot be expected to be 
better than about 10 percent. The exact value of 
the asymptotic constant Ai is not very critical. A 
variation of 10 percent in A z affects the fit by more 
than the expected 3 percent error in fitting only on 
the tail of the distribution. The results in this 
region are not expected to be very precise anyway 
due to neglecting the range straggling. The values 
of the energy dissipation function near the source 
are of low precision because of the very slow con- 
vergence of the Legendre series in this region and 
also of low precision in the backward hemisphere for 
all penetrations because these values result from 
taking differences of nearly equal numbers. 

An alternative double-fitting procedure would be 
to use the spatial moments in cylindrical coordinates 
(2, p). This presents an additional difficulty due to 
the fact that while the fitting in z is straightforward 
the asymptotic form to use in the fitting in p is not 
known. 
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